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Surfaces sputtered by ion beam bombardment have been known to exhibit patterns whose be- 
havior is modeled with stochastic partial differential equations. A widely accepted model is the 
Cuerno-Barabasi model which is robust in its predictions of sputtered surface morphologies. An 
understanding of the factors responsible for such surface topographies can be achieved by using 
scaling arguments on the stochastic model. For such explanations, knowledge of the coefficients is 
crucial. The more so since these vary with different materials, the sputtering process itself gener- 
ates non-equilibrium surfaces within some finite timescale, and the implication of recent results of 
surface topographies unexplained by the continuum theory. We calculate and study these coeffi- 
cients as functions of the sputtering parameters for yet unreported cases of anisotropic ion energy 
distribution within the sputtered material. Consequently, we present phase diagrams for the signif- 
icant case of anisotropic ion straggle. We observe shifts in the phase boundaries when the collision 
cascade geometry rotates, and we also found saturation behavior in the diagrams; in which case the 
boundaries become independent of the penetration depths. Our results indicate a possible origin 
of yet unexplained nanodot topographies arising from oblique incidence ion etching of amorphized 
surfaces. 



I. INTRODUCTION 

Surface sputtering is a process by which materials are 
removed from the surface of a solid through the impact 
of energetic particles. There are many sputtering tech- 
niques, some of which are for industrial uses. It is a 
widely applicable technology with remarkable level of so- 
phistication. Sputtering can be used for surface analysis, 
depth profiling, surface cleaning, micromachining, depo- 
sition, surface coating, semiconductor doping, etching, 
magnetic storage technology, design of nanostructures on 
a surface and many morei - — 

The sputtering process creates patterns on the surface 
at nanometer length scales; hence there is the need to 
study the surface morphology for knowledge and control 
of the pattern of these nanostructures exhibited on the 
surface when it is sputtered. The periodic ripple patterns 
observed for off-normal incidence sputtering were found 
to have either of two orientations, depending on the angle 
of incidence of the impinging energetic ions bombarding 
the surfaceji£~— The ripples are oriented perpendicular to 
the projection of the ion-beam direction onto the surface 
plane for small incidence angles, whereas they are ori- 
ented parallel to this projection for grazing angles. How- 
ever, a number of experiments on surface sputtering have 
observed only a kinetic roughening of the surface, or nan- 
odots, instead of the periodic ripple patterns observed by 
others,! -^ 

Bradley and Harper (BH) were the first to theoreti- 
cally account for the ripple morphology on amorphous 
substrate using the linear continuum theory^ which en- 
tails the theoretical modeling of the surface as a contin- 
uum of points whose evolution with time is governed by a 
deterministic partial differential equation. Their theory 
provides an understanding of ripple formation and depen- 



dence of orientation on ion incidence, but predicts an ex- 
ponential increase of the ripple amplitude. Experiments, 
however, observed rough surfaces at larger lengthscales, 
and a saturation of the surface width. 10 ' 11 

An improvement of the linear BH theory is the widely 
acceptable Cuerno-Barabasi (CB) model which is robust 
in its predictions of sputtered surface morphologies 
CB model is a stochastic nonlinear continuum model 
which considers nonlinearities in addition to the linear 
theory. CB derived this stochastic non-linear equation, 
that mimics the surface evolution with respect to time, 
by working in the laboratory frame of reference, such that 
the coefficients appearing in the equation are functions 
of the actual experimental sputtering parameters. An 
understanding of the factors responsible for such surface 
topographies predicted by CB can be achieved by using 
scaling arguments on the stochastic model based on the 
relative signs of the coefficients; which are the weighting 
factors in the continuum models 

However, experiments by Fackso and co-workers on 
normal incidence sputtering of GaShlS and S i 12 ' 19 ob- 
served nanodots instead of ripple morphologies or ki- 
netic roughening. Moreover, recent results of the dis- 
crete theory in Ref. and [U from Monte Carlo sim- 
ulations of the time evolution of the discretized surface, 
predicts the presence of such nanodots for oblique inci- 
dence sputtering as well. This has been recently con- 
firmed by experiments^ Much is still unclear about the 
continuum theoretical explanations of these recent to- 
pographies. Thus, knowledge of the surface topographic 
weighting factors, of the stochastic partial differential 
equation that describes the surface evolution, and their 
comparative analysis is crucial for an understanding and 
explanation of the predicted surface topographies and 
morphology and for knowledge of how these differ from 
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II. CONTINUUM THEORY 



Energetic ions 




Energy distribution 



For clarity in what follows we provide a schematic di- 
agram, of the ion bombardment of the substrate in a 
sputtering process, below in Fig. [1] 

The linear theory of BBi£ is given as Equation Q] while 
the non-linear theory of Cuerno Barabasiii is given as 
Eq. [2 Equation [5] is an improvement on [T] with the 
introduction of the non-linear terms. 
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FIG. 1. A schematic diagram of the ion sputtering process. 



the continuum theory predictions. Since the continuum 
theory calculations have, however, not yet been per- 
formed for values of the sputtering parameters of the 
order of those used in the simulations it is crucial that 
this be done in order to ascertain whether, among other 
things, the nanodots predicted and observed for oblique 
incidence sputtering constitute a breakdown of the the- 
ory. 

Hence, in this paper we calculate and comparatively 
study these sputter coefficients as functions of the sput- 
tering parameters for the case of anisotropic distribution 
of the energy of the impinging ion within the sputtered 
material. In which case, we have performed calculations 
for values of the order of those in the simulations. But, 
in this first instance, we restrict ourselves mainly to a 
comparison of the relative strengths of the anisotropic 
cascade parameters. Our results indicate shifts and satu- 
ration of phase boundaries due to the relative strengths 
of the perpendicular cascading indicators, and suggests a 
possible origin of yet unexplained nanodot topographies 
arising from oblique incidence ion etching of amorphized 
surfaces. 

The rest of this paper is organized as follows. In the 
next section we discuss the linear and the nonlinear con- 
tinuum theory. In section UTTl we present and discuss our 
results of the phase boundaries for the yet unreported 
anisotropic cases. Finally, we give our conclusions in sec- 
tion HV1 
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h is the surface height at position r and time i, and 
vo is the erosion velocity of a flat surface. (3 X (or v x ) 
is the surface tension coefficient along the x direction, 
which accounts for the instability created by the sput- 
tering process in which surface depressions (troughs) are 
eroded in preference to protrusions (crests). \ x is the 
nonlinear coefficient along the x direction, which arises 
as an improvement of the linear theory to account for 
the stability of the ripple amplitude and other nonlinear 
effects observed in experiments. D is the surface diffu- 
sion coefficient, which takes the thermal hopping of sur- 
face atoms into consideration; such migration of surface 
atoms is the counteracting process tending to restore the 
surface tension instabilities. 

The surface patterns observed arises from the right bal- 
ance of the interplay between the surface tension insta- 
bilities arising from the erosion process and the stabi- 
lizing effects of surface diffusion. r\ is the noise term 
representing the random nature of the sputtering of sur- 
face material; it is taken to be a random variable with 
a Gaussian distribution and zero mean. The difference 
between /3 X and u x is that the former is a function of the 
local sputtering parameters in the local surface frame of 
reference, whereas the later is a function of the experi- 
mental sputtering parameters in the laboratory frame of 
reference. Thus, apart from Equation [5] being a more re- 
alistic theoretical picture of the sputtering process, it is 
also the most relevant in our investigations here for direct 
comparison with the experimental conditions and, hence, 
results. The coefficients occuring in Equation [2] (derived 
in Refs. |T3 and 0) are provided for ease of reference in 
the rest of this paper as follows. 
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where 



vj = ft 2 sin 2 6 1 + a 2 cos 2 i 



exp (—a 2 a 2 cos 2 9/2w) 



F is the ion flux and J is the proportionality constant 
between the power deposition and the rate of erosion. 
For convenience we have set FEJ = 1, which implies the 
coefficients are in units of FEJ. Cuerno and Barabasi 
presented results of phase diagram calculations for the 
isotropic case for simplicity and convenience in their in- 
troduction of the theory*^ They found three different 
scaling regimes defined by the relative signs of the sur- 
face tension and nonlinear coefficients. These scaling 
regimes are applicable to an understanding of the large 
length scale (of the order of tens of m) topographies as 
found in experiments in which topographies have been 
presented for lower resolution of the probe instruments 
[often atomic force microscopes (AFM) or scanning tun- 
neling microscopes (STM)]. 

In this paper we have interpreted our results for the 
anisotropic case in terms of these three scaling regimes 
when nonlinearities are relevant. Recent work in Ref. 
I22I found a much larger number of possibilities for the 
anisotropic case. Following Cuerno and Barabasi, the 
regions to be encountered below are defined by the signs 
of the linear and nonlinear coefficients as follows - I: v x < 
0, v y < 0, < 0, X y < 0; II: v x < 0, v y < 0, \ x > 0, X y < 
0; III: v x > 0,is y < 0, X x > 0,X y < 0. For regions I 
and II, orientations of the ripples formed are along the 
x direction as this presents the highest magnitude of the 
modulus of the surface tension coefficients. Region III is 
characterised by orientations along the y direction since 
in this case v x > 0. 



III. 



RESULTS AND DISCUSSION 



The results presented here are obtained from numerical 
computations of the equations of section [II] for a ranging 
between and 3, and for values of ranging between 
and 90°; where we obtain a set of values for each of 
v x j y and X x / y , and perform a comparative study of their 
relative signs for the phase plots. We have not consid- 
ered higher values of a in order to be able to display the 
saturation behaviour, which we found to depend on a, 
without considering values of the collision cascade pa- 
rameters higher than those in the simulations. 
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FIG. 2. Phase diagram for the isotropic case a = fj, = 7. The 
phase boundaries are found to saturate. The regions are as 
defined and explained in the text. 



Shown below are contour plots of the obtained data 
with which we shall make comparison and then make 
conclusions. First, we consider higher values of a and 
fi that have not yet been reported for the isotropic case 
introduced in Ref. Qj], and found the saturation behavior 
presented in Fig. [2] below for a = fi = 7.0. As can be 
seen from this figure, the phase boundaries become inde- 
pendent of the penetration depth, a, for wider ellipsoid 
of ion straggle within the substrate (see Fig. [1]) . 

This implies that, as values of a and \i increases the 
surface ripple patterns and their orientations, as defined 
by the regions, are observed with the phase boundaries 
depending on 9 and independent of the depth (Fig. [2]) . 
In the rest of the paper we present our results for the 
anisotropic case, a + 1 /1. In Figs. [3] and |H we present the 
large length scale sputtering behaviour, which consists 
of the regions I, II, and III as defined above in section 
2. The short length scale behaviour, discussed below, is 
presented in Figs. [5] and [5] 

A comparison of the phase diagrams of Figs. [5] to @] in- 
dicates that the phase boundaries of region two centered 
around 9 31° in Fig. [2] shift to higher values of 9 as 
the relative size of the longitudinal stragle with respect 
to the lateral straggle parameter /i increase, and to lower 
9 < 31° as the relative size of a decrease, such that the 
dominance of the patterns of regions I extend to higher 
or lower sputtering angles, respectively. With decreasing 
relative size of cr, we found a shift of region I upwards 
due to constraints imposed by the increasing dominance 
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FIG. 3. Phase diagram for the anisotropic case a = 3.5, 
fi = 1.5 [top]; and vice-versa. 



FIG. 4. Phase diagram for the anisotropic case a — 6, fi - 
[top]; and vice- versa (b) 



of the topography defined by region III. Thus, surface 
patterns of region III now dominate almost the entire 
morphology for experimentally significant incidence an- 
gles, and penetration depths. And, topography defined 
by region I, or region II, is almost irrelevant for common 
sputtering angles. 

According to Cuerno and Barabasijil the three regions, 
found and presented thus far, are applicable to the large 
length scale properties of the surface morphology, and in- 
dicative of the three scaling regions characteristic of sput- 
tered surfaces at such length scales. As pertains to ripple 
morphologies alone, therefore, there are only two possible 
regions either one having an orientation that is perpen- 
dicular to the other, as explained in section [XT] Such 
ripple topographic regions account for the short length 
scale behaviour. In the rest of the paper, we restrict our- 
selves to this characterisation of the surface morphology 
in terms of ripple topographies and present results for 
some anisotropic cases. The regions in the phase dia- 
grams presented below are defined as follows, with the 
usual orientation characteristics: I. v x < v y < 0; II. 
v y < v x . 

In Figure O top, we observe that patterns on the sur- 
face are more of the orientation defined by region I than 
that defined by II. For angles less than 6c\ ~ 54°, the rip- 



ple orientation is definitely along the x-direction, and for 
incidence angles greater than 9c 2 ~ 68°, the ripple orien- 
tation is along the y-direction for this case of anisotropic 
collision cascade geometry. For angles within the range 
0c x < < 9c 2 1 the ripple orientation depends on the an- 
gle of incidence and penetration depth of the impinging 
ion within the substrate as in the figure. This interpre- 
tation of the result of Fig. 8 applies to Figs. [5] (bottom) 
and [HI in which cases we also have anisotropic geome- 
tries of the ion straggle. We will just note here that the 
difference is due to the anisotropy, and that changes in 
the boundaries shown in the figures are due to the rela- 
tive magnitudes of the cascade parameters as explained 
above. Note the saturation behavior displayed in the 
phase diagrams of Fig. [51 



IV. 



CONCLUSIONS 



We have performed phase diagram calculations for 
practical anisotropic colision cascade parameters yet un- 
reported in order to gain clearer insight into the view- 
point of the continuum theory as regards the origin and 
theoretical explanation of recent nanodot topographies 
observed for off-normal incidence sputtering. In any case 
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FIG. 5. Short length scale phase diagram for the anisotropic FIG. 6. Short length scale phase diagram for the anisotropic 
case a = 2, /u = 1 [top]; and vice- versa. The regions observed case a — 2, fj, — 6 [top]; and vice- versa, 
are as defined in the text. 



the anisotropic inclosing geometry of the ion straggle 
within the substrate is the most likely in reality but un- 
reported as presented here largely because the isotropic 
cases that have been reported so far have been mainly 
for simplicity and ease of presentation of the theory. In 
this first instance, we have considered the three scaling 
regions defined by Cuerno and Barabasi in their intro- 
duction of the widely acceptable form of the continuum 
theory. 

We found shifts in these phase boundaries as the rela- 
tive magnitude of the collision cascade parameters vary, 
with saturation for higher values of these parameters. 
The saturation behaviour has been found for values of the 
collision cascade parameters higher than the penetration 
depth. Specifically, we observed that when the relative 
size of the longitudinal straggle a is greater than unity, 
the patterns and large length-scale scaling behaviour of 
region I dominates while at the same time the phase 
boundaries straighten from their curved outline until sat- 
uration, with increasing a. On the other hand, when the 
relative size of a is less than unity, patterns (and orien- 
tation) and scaling of region III dominates while at the 
same time the influence of region I reduces with increas- 
ing \i. The saturation behaviour of the phase boundaries 
is also found in this case. 



This exchange of patterns and orientations strongly 
suggests a superposition of several patterns and orienta- 
tion rapidly or gradually (as the case may be) occurring 
on the topmost surface layer downwards. This may be 
the factor responsible for the recently reported nanodot 
topographies predicted for oblique incidence ion sputter- 
ing. Hence, while the results are not yet conclusive about 
nanodot formation, they at least indicate that the nan- 
odot topographies may not constitute a breakdown of the 
continuum theory. 

Through our comparative study of the weighting fac- 
tors in the stochastic partial differential equation that 
describes the surface evolution, we have been able to pro- 
vide these topographic phase diagrams predictive of the 
large length scale and the short length scale behaviour to 
be expected in experimental investigations that will al- 
most always obtain anisotropic collision cascades in their 
surface sputtering. Experimental investigations of the 
scaling and topographic regimes reported here will consti- 
tute an important step in an understanding of the scaling 
exponents characterizing the large length scale regions. 
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